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ABSTRACT 

We perform binary stellar evolutionary calculations following the simultane- 
ous evolution of both stars in the system to study a potential progenitor sys- 
tem for the Type lib supernova 2011dh. Pre-explosion photometry as well as 
light-curve modeling have provided constraints on the physical properties of the 
progenitor system. Here we present a close binary system that is compatible 
with such constraints. The system is formed by stars of solar composition with 
16 Mq + 10 Mq on a circular orbit with an initial period of 125 days. The pri- 
mary star ends its evolution as a yellow supergiant with a mass of ^ A Mq, a final 
hydrogen content of ~ 3 — 5 x 10~^ Mq and with an effective temperature and 
luminosity in agreement with the HST pre-explosion observations of SN 201 Idh. 
These results are nearly insensitive to the adopted accretion efficiency factor f3. 
At the time of explosion, the companion star has an effective temperature of 22 
to 40 thousand Kelvin, depending on the value of /3, and lies near the zero age 
main sequence. Considering the uncertainties in the HST pre-SN photometry the 
secondary star is only marginally detectable in the bluest observed band. Close 
binary systems, as opposed to single stars, provide a natural frame to explain 
the properties of SN 2011dh without the need to fine tune any parameter. 

Subject headings: stars: evolution — binaries: close — supernovae: general — 
supernovae: individual (SN 2011dh) 
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INTRODUCTION 



Core-collapse supernovae (CCSNe) are the explosive end of massive stars with 
^ZAMS ^ 8Mq. There is a diversity in the spectroscopic and photometric properties of 
CCSNe which are mainly related to the ability of the progenitor to retain its outermost 
layers. Type II SNe, with clear H lines in their spectra, represents the case where a thick 
H envelope is kept before the explosion. Type lb SNe, with no H lines but with clear He 
lines, have lost their H envelope but not the He layers. Finally Type Ic SNe, with no 
H and He lines in the spectra, represent a more extreme case where not only the the H 
but also the He envelopes are likely lost before the explosion. There are also transitional 
objects between these different types. One example is that of Type lib SNe, which show H 



lines at ear 



Filippenko 



y tim es but then the spectrum is transformed into that of typical SNe lb (see 
19971 for a classi fication scheme). Type lib, lb and Ic objects are collectively 



called striped-envelope SNe (IClocchiatti et al. 



19961 ). 



Progenitor models of SNe lib comprising a helium star surrounded by a very thin 
hydrogen- rich envelope (of < 1 Mc^) have been successful to explain the observed light curv es 



(LC) 



and spectral features 



Shigeyama et al.l 11994) : iWoosley et al. 



1994 



Blinnikov et al. 



19981 ). However, it is not clear which is the mechanism responsible for the removal of the 
outer envelope before the explosion. One possibility is strong winds that occur in massive 
stars with Mzams ^ 25Mq. Alternatively, in close binary systems (CBS) stars are expected 
to exchange mass providing an efficient mechanism to allow for the removal of outer layers. 



2008; 


Claevs et al. 


2011; 


Smith et al. 


2011) 



^OGB is member of the Carrera del Investigador Cientifico de la Comision de Investiga- 
ciones Cientificas de la Provincia de Buenos Aires (CIC), Argentina. 
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Additional support for the binary sce nario in SNe 



companion for the famous SN lib 1993 J (Maund et ah 



by pre-explosion photometry (lAldering et al 



1994 



lb comes from the detection of a hot 



20041 ). This was initially suggested 



)■ The LC and evolutionary models of 



SN 



1993 



993J were also in favo r of the binary channel (INomoto et al. 



Wooslev et al 



I993I : 



Podsiadlowski et al. 



19941). Some evidence for a companion was also reported for another 



SN lib 2001ig f,Rvder et al. 



20061 ). 



The Type lib SN 201 Idh was recently discovered in the nearby galaxy M51 attracting 
the attention of many observers bec ause of its proximit y and brightness. It was discovered 



almost im mediately after 



archive 



emissio n (jSoderberg et al. 



Maund et al. 



explos ion (lArcavi et al. 



201 ll ). It showed early radio and X-ray 



20121 ) ■ Using pre-exp 



(120 111 ) and 



Van Dyk et al. 



osion images obtained from the HST 



( I2OIII ) detected a source at the location of 



SN 2011dh. They derived similar values of luminosity and effective temperature for pre-SN 
source. The object was consistent with a yellow supergiant (YSG) star with a radius R 
~ 270 Rq and without any clear evidence of a companion star contributing to the observed 
spectral energy distribution (SED). 

At present there is a controversy in the literature as to whether the YSG is the actual 
progenitor of S N 201 Idh. Some authors have suggested that the exploding sta r should be 



more compact (lArcavi et al 



2011 



Soderberg et al 



2012 



Van Dyk et al 



2OIII) based on 



(1) a simple comparison between the early light curve (LC) of SN 201 Idh and SN 1993 J; 

(2) a discrepancy between the temperature derived from an early-time spectrum and that 
predicted by an analytic expressions for an extended progenitor; and (3) the large shock 
velocity derived from radio observations. 

Recently, we have performed a d etailed hydrodynam ical modeling of SN 2011dh 



using stellar evolutionary progenitors (iBersten et al 



2OI2I ). These models indicate that 



observations are compatible with a helium star progenitor of a mass near 4 Mq surrounded 
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by a thin hydrogen-rich envelope (~ 0.1 Mq) with a radius of ~ 200 i?© (similar to that of 
the detected YSG star) that underwent an explosion with an energy of 8 x 10^° erg that 
synthesized 0.063 Mq of ^^Ni. Such large radius values are needed to reproduce the early 
light curve of SN 201 Idh without contradicting the temperatures derived from the spectra. 
In addition, our hydrodynamical modeling rules out progenitors with He core masses larger 
than 8 Mq, which corresponds to Mzams ^ 25Mq. 

It is very difficult for a single star to reach these pre-SN conditions. The existence of 



a strong win d capable of 



Mq or more ( iHeger et al 



removing most of t 



2003 



le env elope requires a massive star of ~ 25 



Georgy et al. 



20091 ). which is in contradiction with the LC 



models. Moreover, in order to retain a thin hydrogen-rich layer, the mass loss rate would 
have to be on a very narrow interval. These facts strongly suggests that the progenitor of 
SN 201 Idh should be a component of a binary system. 



However, a recent work by 



Georgyl ( 120121 ) proposed that single YSG stars such as the 



one detected at the location of SN 201 Idh are plausible SN progenitors. This is based on 
stellar evolution calculations of stars with main sequence masses of 12-15 Mq under the 
assumption of an increased mass-loss rate several times above the standard values. However, 
no physical explanation is given for such an increas e. Moreover, recent observations 



suggested that the st andard stellar wind mass-loss fjde Jaeer et al 



considerably reduced ( iBouret et al 



2005 



FuUerton et al. 



19881) needed to be 



20061 ). Although eruptive mass 



loss duri ng late stages of stellar evolution may be an alternative way to lose larger amounts 



of mass ( ISmith &: Owocki 



20061 ). some observational evidence in the spectra of the SN due 
to the interaction between the ejecta and previously expelled GSM is expected as in the 
case of SNe Iln. This was not observed for SN 2011dh. 

The aim of this work is to show the plausibility that the progenitor of SN 201 Idh 
was part of a close binary system (GBS) with properties compatible with the pre-SN 
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observations and the results of LC modeling. The observational properties of the remaining 
companion star are discussed in anticipation of future detections. Although we do not 
perform a complete exploration of the parameter space (stellar masses, initial orbital period, 
and mass-transfer efficiency we show that our results are robust if we consider moderate 
changes of the initial conditions. 

The remainder of this paper is organized as follows. In Section |2] we present a brief 
description of our binary stellar evolution code paying special attention to the characteristics 
that enabled us to compute pre-SN models. In Section [3] we present the main results of 
this paper regarding the adopted binary configuration (§ 13. ip . evolutionary calculations 
(§ 13.21) as well as the spectra of the components at the moment of the explosion (§ 13. 3p . In 
Section H] we present a discussion of our results and finally, in Section [5] we provide some 
concluding remarks. 



2. THE STELLAR CODE 



In order to compute th e evolution of the CBSs quot ed above we shall employ a code 



similar to that described in 



Benvenuto &: De Vitd (120031), adapted for the case of massive 



stars. Briefiy, this is a Henyey code that when the star reaches Roche Lobe Over Flow 
(hereafter RLOF) conditions it solves implicitly not only the whole structure of the donor 
star but also the Mass Transfer Rate (hereafter MTR) Mi in a simultaneous, fully implicit 
way. Such a procedure has been found to largely improve numerical stab ility as compared 



to algorithms that compute the MTR explicitly (Biining fc Ritter 



20061), enabling us to 



reach pre-SN condition meanwhile the donor star is losing mass at an appreciable rate with 
a detailed, continuous and convergent sequence of stellar models. In detached conditions, 
the code works as a standard Henyey scheme. 
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In order to adapt the code to the case of massive stars we have 



nuclear reactions from the compilation of 



Caughlan fc Fowlerl (119881 ) and rewritten 



incorporated several 



the difference equations to largely improve numerical stabilit y. Also, we i n corpo rated 



semiconvection following the diff usive approach 



first applied to massive stars by 



Langer et al. 



prese nted in 



Langer et al. 



(119831) and 



(119851 ). This formulation of diffusive 



se miconvection depe nds on an efficiency parameter asc for which we assumed a5C7 = 1 as 



m 



Yoon et al. 



(|2010[ ). In this paper we shall ignore the effects due to overshooting and 



rotation. 



In semidetached conditions we compute Mi following the prescription given by 



Ritter 



mm 



Ml = -Mo exp 



Ri — Rl 
H„ 



(1) 



where Ri is the radius of the donor star; R^, the r adius of a sphere with a volume equal 



EggletonI (11983[ ): Mq > is a smooth 



to the corresponding Roche lobe is computed following 
function of Mi and M2 whereas Hp = —dr/dln{p) (where r and p are the radial coordinate 
an d pressure res pectively) is the photospheric pressure scale height. For fu rther details 



see 



de Jager et al 



Ritterl (119881 ). In detached conditions we consider stellar wind following 
(jl988l ). As we shall discuss in § [3l we consider different accretion efficiencies, /3, of the 
material transf erred by the donor to the sec ondary, i.e. M2 = —(3 Mi, and evolve the orbit 



as described in 



Benvenuto fc De Vitol toom . 



After computing the evolution of the donor star, we evolve the com panion star takin 



Neo et al. 



( 11977h 



into account the accretion rate it receives from the donor as described in 
We do not consider mass loss from the accreting star. As we ignore any effect of the 
secondary star on the donor one other than imposing a limit on its volume, this is an usual 
assumption in this type of problems. 
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It is well known that in CBSs like those we are interested in here, the material accreted 
by the secondary star may have a mean molecular weight higher that the corresponding 
to its outer layers. T his may lead to an unsta ble situation that induces the so - called 



thermohaline mixing (IKippenhahn et al. 



19801). In the present version of our code we 



have not included th e rmoh aline mixing yet. In any case, it is worth to remark that 



Stancliffe &: Eldridgd (120091 ). have shown that ignoring thermohaline mixing has a minor 
effect on the evolutionary track of secondary stars of CBSs similar to those studied here. 
Therefore, we expect that the main conclusions of this work not to be affected by neglecting 



this phenomenon. 



3. THE BINARY MODELS FOR SN 2011dh 
3.1. The initial configuration 



Van Dyk et al. 



of the SN 201 


Idh ( 


Maund et al. 


2011 


(Bersten et al. 


2012 


) impose strong 



constraints about our election of the binary parameters, mass of the primary, of the 
secondary and the initial period. Here we briefly discuss the motivations of our election of 
these parameters before we present our results. 

If the object detected in the pre-explosion image is indeed the progenitor of SN 201 Idh 



and it belongs to a binary system, th e primary (donor) s tar n eeds to have a 



compatible with the value derived by 



Maund et al 



( 1201 ih and 



Van Dyk et al 



uminosity (L) 



(120 111 ) for the 



progenitor candidate, i.e. log L/Lq = 4.92 ± 0.2. This luminosity is an indication of the He 
core mass at the moment of the explosion wh i ch in turn is related with the mass of the star 



on the main sequence (Mzams) 



Maund et al 



( 1201 ll ) derived Mzams = 13 ± 3 Mq comparing 



the luminosity with the end points of the evolutionary tracks of single stars. Note that they 
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did not use the color information of the progenitor to derive this mass bec ause uncertaint ies 
are ex pected in the color due to unknown mass-loss history. Meanwhile 



Van Dyk et al.l 



( 120111 ) derived Mzams = 17 — 19 using both L and the effective temperature (Tcs) 
derived from colors and choosing the track that best matched these values, although this 
point does not corresp ond with the f i nal p osition of the single star at the end of the 



evolution. In addition, IBersten et al. 



(120121 ) derived a helium core of ~ 4Mq from the light 
curve modeling of SN 2011dh and firmly ruled out progenitors with Mzams > 25 Mq. Here 
adopt an initial mass of 16 Mq for the donor star that is well within the ranges provided in 
previous studies. 

After assuming a mass value for the donor star, we still have to consider the mass of 
the secondary star as well the orbital period of the binary. Our choice should be guided by 
the fact that pre-explosion observations indicate that the observed portion of the spectrum 
is compatible with a single source. Obviously, the secondary has to be less massive than the 
donor star but we have to distinguish between two cases: i) the mass ratio is close to one 
and ii) the mass ratio appreciably differs from one. 

Let us consider the case i): If the masses of the stars differ in (say) few tents of solar 
mass, the secondary star would be able to exhaust core hydrogen before explosion. This 
object would receive material coming from the donor star when it is on the Hertzsprung 
gap where the shell nuclear burning around the core takes place. Calculations available in 
the literature as well as our own test (see § H]) indicate that, after accretion, such object 
appears as an overluminous B supergiant with a T^g in between that of the observed YSG 
and that of the ZAMS (see, e.g.. Fig. 5 of Claeys, et al. 2011 and our Fig. [15]). In this 
case some evidence of the secondary should have been detected in the HST pre-explosion 
photometry as in the case of SN 1993J. 

Case ii) results more natural. If the secondary star has a mass appreciably lower than 
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that of the donor, it will be still undergoing core hydrogen burning at the moment of 
the explosion of the primary star. Therefore, at the moment of the explosion and in the 
HRD, we expect the secondary to be close to the ZAMS. The object will remain hot and 
will emit most of the flux in the UV. As it is well known the more massive the object the 
greater its luminosity on the ZAMS and therefore the greater its effect on the pre-explosion 
photometry. Thus, the mass value chosen for the secondary star should be low enough to 
have remained almost undetected. Note, however, that if the mass of the secondary is much 
lower than 10 Mq, the Kelvin-Helmholtz timescale increases appreciably and it is very 
likely that the system reaches common envelope conditions. Therefore, we adopt a mass of 
10 Mq for the secondary. 

Even with a mass of 10 Mq if the secondary object were able to undergo a conservative 
mass transfer and the primary star ends its evolution with a mass of w 4 Mq as required by 
LC modehng, the secondary would end its evolution with fa 22 Mq on the ZAMS. This would 
produce a very bright object that may not be compatible with the pre-SN observations. 
One possibility that we will explore in the next section is that the secondary captures a 
fraction /3 of the material transferred, therefore its final mass will be of ~ 10 + 12 /3 Mq. 
Note that the accretion efficiency, (3, is one of the most uncertain parameters in binary 
stellar evolution. Another parameter usually employed in the treatment of orbital evolution 
in binaries is the specific angular momentum a of the material lost by the system in units of 
the angular momentum of the primary star. We shall assume a — 1 throughout this paper. 

Regarding the initial period, this has been evaluated in order for the pre-SN donor star 
to fall inside the error box given by the L and T^s estimated for the progenitor candidate 
of SN 2011dh. Suppose that the donor star is losing mass by RLOF at the moment of 
the explosion, i.e. that the the size of the Roche Lobe is approximately equal to the 
radius of the donor star. This radius can be determined using log L/Lq — 4.92 ± 0.2 and 
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logTcff = 3.78 ± 0.02 flMaund et al 



2011 



Van Dyk et al. 



20111 ) leading to values of ^ 270 



Rq. If the initial masses are 16 + 10 Mq and the final masses are 4+12 Mq, the the final 
orbital semiaxis is ~ 900 Rq and the final period at the moment of the explosion can be 
estimated to be ^ 8 day s. Finally, the initial period can be calculated using Eqs. 6-8 of 



Podsiadlowski et al. 



( I2OO2I ) which relates the initial and final orbital semiaxes as a function 
of the initial and final masses of the system. This leads to initial periods of ~ 120 days. 
Note that if the luminosity is due to the internal structure of the star, then at pre-SN 
conditions a given orbital period will correspond to a given effective temperature: the larger 
the period the lower the effective temperature. In the next section we analyze in detail a 
system with an initial period of 125 days consistent with this first approximation. Other 
values are discussed in section [51 



3.2. Evolutionary Results 

In the previous section we discussed our election of the binary parameters adopted 
to study a possible progenitor for SN 2011dh. Here we present our results for a CBS of 
solar-composition stars with masses of 16 Mq + 10 Mq on a circular orbit with an initial 
period of 125 days. As stated in § 13. Ij we analyzed different values of the mass-transfer 
efficiency, (3 = 0.00,0.25,0.50,0.75, and 1.00. We computed the evolution of both stars 
starting with ZAMS models up to core oxygen exhaustion. The main results of these 
calculations are presented in Figures [H to [121 

As an extreme case, we show in Fig. [H the evolution of the donor and companion 
stars assuming /3 = 0.00 (fully non - conservative mass transfer). This CBS undergoes 
class B mass transfer: the donor star fills its Roche lobe with a mass of 15.54 Mq well after 
core helium ignition (the central abundance of helium at the onset of the first RLOF is 
^He = 0.32). With respect to the stellar mass, this is the main RLOF because the donor 
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star detaches from its lobe when has only 4.54 Mq. Remarkably, the RLOF takes only 
7.6 X 10^ y, which implies a mean MTR (M) of (M) = 1.44 x 10"^ M^/y. However the 
maximum mass transfer rate M^ax is of M^ax = 1-90 x 10^'^ ^q/u (see below for further 
details on MTRs). At detaching from its Roche lobe, the donor star is still undergoing 
core helium burning and has increased its central abundance due to semiconvective mixing 
(X//e = 0.39). The physical agent that sets the end of the RLOF is the outer layers 
hydrogen abundance Xh\s that fallen down from its initial value of Xh\s = 0.70 to 
Xh\s = 0.49. Thus the stellar envelope is no longer able to support its very large radius 
(192 -Rq), and starts a fast contraction performing a blue-wards loop on the HR diagram. 
During this loop core helium is exhausted. Soon, due to semiconvective mixing and nuclear 
shell burning the star undergoes a second and very brief RLOF followed by another (smaller 
as compared to the previous) blue-wards loop. Finally, the star swells and again goes 
through a RLOF up to its final explosion. During this RLOF the star undergoes carbon, 
neon, oxygen and silicon core burning (we did not computed the silicon burning stage). 
During this final RLOF the mass transferred from the donor is a relatively small amount of 
0.27 Mq on 3.57 x 10^ y ((Af ) = 7.95 x 10"^ Mq/v and Mrnax = 2.92 x 10"^ M^/y). Notice 
that these MTRs are far lower than the corresponding to the initial RLOF. 

The evolution of the central point of the donor star in the temperature-density plane 
is presented in Fig |2j There, it can be clearly noticed the nuclear activity present in the 
central part of the star during each RLOF. It is remarkable that during the first RLOF 
central density and temperature remain almost constant. 

In this case of /9 = 0, the secondary star does not accrete any material and evolves 
as it were an isolated object. As the main sequence lifetime of a 10 Mq star is far longer 
than that corresponding to the donor star, the companion star suffers from a very small 
excursion in the HRD up the moment of explosion (see Fig. [1]) 
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The evolution of the components of the CBS for the cases of /3 = 0.25, 0.50, 0.75 
and 1.00 is shown in Figs. ISE respectively. An inspection of these figures is sufficient to 
realize that the evolution of the donor star is almost independent of the value of /3 while 
the secondary is strongly dependent on it. In Fig. [7|we show the evolutionary tracks of 
the donor star corresponding to the five values of /3 considered in this paper. They are 



remar kably similar. This be 



lavior resembles the results found in the case of low mass 



CBSs flDe Vito fc Benvenuto 



2OI2I) 



In Fig. [H] we show the evolutionary tracks of the companion, accreting star as a 
function of /3. There, for comparison, we also show the ZAMS corresponding to the initial 
composition of these stars. For the cases of /3 > the final position in the HRD for the 
accreting star is somewhat hotter and overluminous than objects of the same mass on the 
ZAMS. This is partially due to the fact that we have neglected thermohaline mixing in our 
calculations. Note that in all considered cases, the companion star does not fill its Roche 
lobe and no contact configuration is found. Thus, the CBSs studied here do not undergo 
any common envelope episode. 

The evolution of the central point of the accreting star is shown in Fig. [H Notice that 
the excursion of these objects in the density - temperature plane is by far smaller than 
the one corresponding to the donor stars. Accreting stars are not able to exhaust central 
hydrogen in the time spent by the donor star to reach pre-SN conditions 

The evolution of the MTR from the donor star is presented in Fig. [101 The only 
difference in the MTR evolution as a function of /3 is in the time spent by the star from 
the end of the first RLOF to the onset of the second one which is minimum for /3 = 0.50. 
All the cases predict a third mass-transfer episode that lasts until the end of the evolution 
and whose MTR differs markedly from constant. Therefore any inference on the mass of 
the envelope of the donor star at the time of the explosion should take into account the 



- 14 - 



appropriate mass loss in this phase. 

From the behavior of the evolution of the MTR it is also posible to understand 
quantitatively the evolution of the companion star. During the first RLOF the MTR is 
large enough (M > lO'^M©/?/) to force the secondary star to swell appreciably, reaching 
low effective temperatures if this material is efficiently accreted (see, e.g.. Fig. Elfor /3 = 1). 
For the range of initial periods considered (~ 125 d) this occurs during core helium burning. 
The first RLOF ends before helium core exhaustion and later the star undergoes a final 
RLOF on which it evolves up to pre-SN conditions. Remarkably, while the MTR in the 
first RLOF may be large enough to inflate the star, calculations indicate that during the 
second and thrid RLOF this is not the case. In fact, in the last RLOF, M is about two 
orders of magnitude lower. So, after the end of the first RLOF and even if the star reached 
very low effective temperature values, it contracts very fast. In the HRD, the secondary 
star moves blue-wards near the position corresponding to an object of the same mass close 
to the ZAMS and remains there until the explosion on the donor star. 



The evolution of the orbital period is n ot very sensitive to t 



seen in Fig. [TTl This is not surprising, since 



le value of /3, as can be 



Podsiadlowski et al. 



(120021 ) showed that the 



evolution of the orbital period of CBS depends only slightly on the accretion efficiency. 
Their analytic prediction for the evolution of the period as a function of the initial and final 
masses and the /3 parameter fit almost perfectly with our numerical results. 

An important result of our calculations is that despite of the fact that the system is in 
semi - detached conditions at the moment of explosion, the donor star retains an appreciable 
amount of hydrogen in its outermost layers as it is shown in Fig. [121 The total hydrogen 
content is again almost independent of the value of /3 and enou gh to account for the H 



lines observed in the spectra of SN 2011dh (jPessart et al 



20111). While a full exploration 



of the total hydrogen content of donor stars at pre-SN conditions for CBSs in general is 
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beyond the scope of the present paper, this result strongly indicates that the total hydrogen 
content should be a function of the initial orbital period: the larger the period the larger 
the hydrogen content. In this sense, the progenitor of SN 2011dh may be considered as a 
transition object. 

For completeness, we present in Table [T] the main characteristics of the pre-SN object 
and in Table |2] those properties corresponding to the accreting, secondary star. 

In summary, we have shown that a system with 16 Mq + 10 Mq and an initial period 
of 125 days, independently of the adopted value of /3, predicts that the primary star ends 
its evolution w ithin the region of the H-R diagram comp atible with the pre-SN photometry 



of SN 2011dh flMaund et al. 



2011 



Van Dyk et al. 



20111 ). Furthermore, at the end of the 



evolution the primary star has a mass of ~ 4 M^^ and a hydroge n content of 3 — 4 x 10 ^Mg 



which is consistent with the LC modeling (IBersten et al 
of SN 201 Idh. 



2012) and the SN lib clasification 



3.3. Pre-Supernova Spectral Energy Distribution 



Our model calculations predict that the primary star ends its evolution with properties 



(L and Teg) compatible with those inferred for 



SN position (IMaund et al. 



2011 



Van Dyk et al. 



:he p re-explosion object located at the 



20111 1. At the same time, the secondary 



star has a luminosity that, depending on the value of /3, could be comparable with that of 
the primary star and therefore produce a detectable effect on the pre-SN spectral energy 
distribution (SED). In all cases the secondary star is found to be significantly hotter than 
the donor star. Indeed, for high enough effective temperatures, the effect of the secondary 
would only be appreciable in the bluest available photometric band. 



We thus study here the effect of the secondary on the SED of the system for different 
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Table 1. Some Properties of the Donor Star at the Moment of Explosion as a Function of 

/3. 





M 


logioCTeff) 


logio(^) 


R 


Mh 


Age 




[Me] 


[K] 


[^0] 


[Re] 


[10-3 Mq] 


[Myr] 


0.00 


4.034 


3.788 


4.886 


245.54 


3.869 


12.66 


0.25 


4.118 


3.790 


4.907 


249.20 


4.441 


12.44 


0.50 


4.118 


3.791 


4.907 


247.73 


4.537 


12.44 


0.75 


4.073 


3.791 


4.889 


242.18 


4.926 


12.53 


1.00 


4.014 


3.786 


4.871 


242.74 


3.465 


12.64 



Note. — Form left to right we tabulate the value of l3; the 
mass; effective temperature; luminosity; radius, total amount of 
hydrogen; and the age of the pre-SN. 
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Table 2. Some Properties of the Secondary Star at the Moment of Explosion as a 

Function of (3. 





M 


logio(Teff) 


logio(^) 


logio(Pc) 


logio(Te) 






[Me] 


[K] 


[^0] 


[g cm~^] 


[K] 




0.00 


10.000 


4.351 


3.810 


0.929 


7.497 


0.426 


0.25 


12.844 


4.473 


4.130 


0.810 


7.513 


0.543 


0.50 


15.688 


4.526 


4.377 


0.729 


7.522 


0.584 


0.75 


18.580 


4.562 


4.581 


0.662 


7.539 


0.597 


1.00 


21.515 


4.592 


4.753 


0.606 


7.548 


0.613 



Note. — Form left to right we tabulate the value of l3; the mass; 
effective temperature; luminosity; and the central values of density, 
temperature, and hydrogen abundance at the moment of explosion. 
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values of /3, and compare this with the HST pre-explosion photometry. To calculate the 
SED for each star of the binary system a s well as the co mposed SED we used atmospheric 
models for solar composition provided by iKuruczl (jl993[ ). The model spectrum of each star 



was obtained by linearly interpolating at the values of Tes and surface gravity (g) given in 
Tables [1] and [21 The observed fluxes were computed b y multiplying the mode ls by {R/dY 



where d is the distance to M51 assumed to be 7.1 Mpc (jTakats fc Vinko 



20061 ) . The sum o 

the SEDs was used to compute synthetic photometry through the HST transmission filter 

Fig. [13] show the resulting SED of each star and their sum for /3 = and 1. Mean 
synthetic and observed fluxes are compa red for all bands. Ob served fluxes were obtained 
from the tabulated magnitudes given by 



Van Dyk et al. 



(120 



if) and zero p oints in the Vega 



Bohlin 



mm . Note that the 



system calculated using the Alpha Lyrae SED provided by 
contribution of the secondary star is significant only for the bluest observed band, F336W. 
The increase in the F336W flux due to the presence of the secondary is of 16%, 22%, 32%, 
45%, and 62% for the cases of /3 = 0.00, 0.25, 0.50, 0.75, and 1.00 respectively. 

For the rest of the observed HST bands the contribution of the secondary is < 6%, 
i.e. for wavelengths A > 4000A the spectra is completely dominated by the light coming 
from the donor, pre-SN star. In this wavelength range the available observations on four 
photometric bands pose no constraint on the properties of the accreting star. On the 
contrary, the bluest filter, with a maximum transmittance at A ~ 3400 A partially detects 
the red tail of the spectrum of the accreting star. 

For the filters with transmittance at A > 4000A the agreement between calculations 
and observations is very good, especially considering that we have not adjusted any 
parameter. For the case of the bluest filter, F336W, the contribution of the secondary to 



^http : //www. stsci . edu/hst/wf c3/ ins_perf ormance/f liters/ 
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the total flux is not negligible. If we consider the measurement uncertainty in this band we 
find that the secondary is detectable only at the 2-o" level in the most extreme case (/3 = 1). 
The contribution is further decreased to the 0.6-a level when completely non-conservative 
mass accretion (/3 = 0) is considered. 

Remarkably, irrespective of the value of /3, at the moment of the explosion, the 
secondary star is so hot that its light has been barely detected. Thus, quite unfortunately, 
the characteristics of the secondary star are very poorly constrained by the presently 
available data. 

4. DISCUSSION 

In the previous section we have shown that a properly chosen binary configuration can 
explain very well the proposed YSG progenitor of SN 2011dh. Our models also predict 
that the total and hydrogen masses of the donor star at the moment of the explosion are 
consistent with the results of LC modeling and the SN lib classification. 

A complete exploration of the parameter space of initial stellar masses and orbital 
periods is not the scope of this paper. Nevertheless, we briefly test that our results are 
robust if we consider moderate changes of the initial parameters. For example. Fig. UM 
shows the sensitivity of the evolutionary track of the primary star on small variations of 
the initial period. Pi between 100 and 200 days, for the same configuration presented in 
the previous section, i.e. 16 + 10 Mq, and f3 = 0.5. Within such range of initial periods 
the donor star ends its evolution as a YSG close to the region of the HRD allowed by 
the pre-SN photometry. The evolutionary track of the secondary is not shown in Fig. [H] 
because it is not sensitive to the initial period in the range under study. 

The effect of changes in the mass of the secondary star on the binary evolution is shown 
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in Fig. [m Three values were considered, M2 = 15, 15.25 and 15.50 Mq. In all the cases, 
we assumed a value of /3 = 0.5 and a mass of the primary star of 16 Mq as in our previous 
models. The initial period was modified to make the end-point of the donor star match the 



region allowed by the pre-SN photometry. Following the criteria described in section 13. 
the adopted initial period for this test was of 50 days. From the figure, it is clear that only 
for the case of M2 = 15.5 Mq, where the mass ratio is closest to one, the secondary star 
moves appreciably away from the ZAMS at the moment of the explosion of the primary. 
This is because the secondary has exhausted its hydrogen core before the beginning of the 
first RLOF. Consequently, the secondary ends with a redder color than that of lower-mass 
stars and it is expected to contribute significantly to the observed pre-SN photometry. 
Indeed, the flux of the secondary in the bluest observed band, F336W, is 2.5 times larger 
than that of the primary star, which would lead to a 8.5 cr detection. For the rest of the 
HST bands the contribution of the secondary is within 6-30%. While not compatible with 
the case of SN 201 Idh, this configuration may be applicable to SN 1993J wh ose pre-SN 



observations showed evidence of a companion star (jStancliffe fc Eldridge 



20091). 



Secondary stars of slightly lower initial masses (15 and 15.25 Mq) remain near the 
ZAMS and end as hot stars of similar luminosity to the primary. Stars of masses between 
10 and 15 will also remain as blue objects near the ZAMS. Because of the high effective 
temperature most of the flux from the secondary is emitted in the UV, away from the 
bandpasses of the pre-explosion imaging. However, if the luminosity is high enough, the 
secondary may produce a detectable effect in the optical range. For instance, with masses 
of 15 and 15.25 Mq, the flux of the secondary in the F336W band would be comparable to 
that of the primary star. The level of detection of the secondary in these cases would be 
close to 3.5 cr. The final luminosity of the secondary also depends on the assumed value of 
p. Lower initial masses could result in similarly high luminosities with values of /3 close 
to unity. Unfortunately, it is not possible to break such degeneracy based solely on the 
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available pre-explosion photometry. Future observations of the possible secondary star after 
the SN fades from sight can shed light on this matter. 

The initial mass of 16 Mq for the primary star leads to a progenitor with the right 
characteristics as derived from pre-explosion photometry, LC modeling and spectral 
classification. Therefore, our choice of initial masses and orbital period is by no means 
unique but it allows to prove that a CBS is a plausible progenitor for SN 201 Idh and other 
SNe lib. Moreover, the right progenitor properties are achieved in a self-consistent manner 
and independently of the detailed initial conditions. 

In the past few years, some progenitors of SN e II have been associat ed with YSG stars , 



e.g. Type IIP SN 200 8cn flEhas-Rosa et al. 



Elias-Rosa et al. 



20091), Type IIL SN 2009kr flFraser et al. 



2010 



20101), and the SN studied here. The explosion of a YSG is not compatible 
with the theoretical prediction of single stellar evolution. Motivate d by thi s appa rent 



Georgvl (120 121 ) studied 



discrepancy and the lack of evidence of a companion for SN 201 Idh, 
the effect of an increased mass-loss rate on the final properties of stars with initial masses 
of 12-15 Mq. Assuming rates several times higher than the standard values, they found 
that it was possible to explain the explosion of single stars of relatively low mass in the 
yellow area of the HRD. In particular, their proposed scenario could explain the progenitor 
of SN 201 Idh. However, no physical explanation was given for such an increased mass-loss 
process. 

Note that the evolutionary tracks of isolated stars are strongly dependent upon the 
details of mass loss. In this sense, explaining the position of the progenitor of SN 2011dh 
in the HRD would require some degree of fine tuning. Moreover, this fine tuning has 
to simultaneously account for the amount of hydrogen necessary to explain the spectral 
evolution. The binary scenario instead provides a self-consistent picture that naturally 
explains why the star remains as YSG for long periods of time and until the explosion. 
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In the case of CBS evolution we find that, irrespective of the mass transfer efficiency, the 
donor star undergoes a final RLOF before igniting carbon and so it is still transferring 
mass to the companion at the moment of explosion. During the final RLOF, the donor 
star transfers a small amount of material. So, the orbit and its Roche lobe enlarge very 
little meanwhile the donor star tends to swell as a consequence of nuclear shell burning. 
This precludes the donor star to attain lower effective temperatures and makes it evolve to 
higher luminosities in a way equivalent to the evolution of giant stars with extended outer 
convective zones. In any case, the observed effective temperature of the proposed progenitor 
of SN 201 Idh strongly indicates the initial orbital period, but no other parameter has to be 
adjusted for the pre-SN object to fall inside the observational error box. The YSG nature 
of the progenitor is thus a direct consequence of close binary evolution. 



Recently IChevalier fc Soderberg) ( l2010l ) suggested a division of SNe lib into compact 
(cllb; R ~ 10^^ cm) and extended (ellb; R ~ 10^^ cm) subtypes essentially based on radio 
LC properties. Compact ob j ects w ere proposed to have smaller hydrogen masses, roughly 



below 0.1 Mq. 



Claeys et al. 



( I2OIII ) used this criterion and their own CBS evolutionary code 
to analyze the range of periods (for Pi > 1000 days) and initial masses needed to produce 
extended SNe lib. However, our calculations show that it is possible to have an extended 
progenitor {R ^ 250 -R0) with a H mass of < 0.1 Mq. Therefore, if the cllb and ellb 
subtypes correspond to physically distin ct progenitors, the division criterion may need to be 



revised. In addition. 



YoonetaL 



(|2010[ ) also analyzed the space of parameters in CBSs to 
produce Type Ib/Ic SNe using very different initial periods of < 8 days, as compared with 
the ones adopted here. They found that some of their models predict a thin hydrogen layer 
of < 0.01 Mq with a compact structure. This channel of production of SNe lib is different 
from the one we have presented and leads to explosions far away from the YSG regime. 
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5. CONCLUSIONS 

With the aim of providing a description of the progenitor of SN 2011dh, we have 
studied the evolution of close binary systems of solar composition stars with masses of 
16 Mq + 10 Mq. We considered an initial period of 125 days and different efficiencies 
of the mass transfer process. We followed the simultaneous evolution of the donor and 
accreting stars from the zero age main sequence up to the oxygen core exhaustion of the 
donor. We found that the donor star, independently of /3, ends its evolution with effective 
temperature and luminosity consistent with the YSG object detected in the HST pre-SN 
photometry. The exploding star has a mass M ^ 4M0, a radius R ~ 250 i?© and an 
outermost layer containing 3 — 5 x 10~^ Mq of hydrogen. This is generally consistent wit 



the ty pe lib classification and the results of LC modeling of SN 2011dh by 



Bersten et al. 



(120121 ). These results are a natural consequence of the close binary evolution and require 



no external adjustment of any physical condition. 

Regarding the accretion efficiency, /3, we found that (1) the evolution of the donor star 
is almost independent of /3 while the secondary strongly depends on it, and (2) the evolution 
of the orbital period, the MTR and the total hydrogen content are almost independent of 
the value of f3. 

Our calculations indicate that the donor star is lossing mass at the moment of the 
explosion with rates that differ markedly from constant. Inferences on the mass of the 
donor star at the time of the explosion should take into account the appropriate mass loss 
in this phase. We also found some indication that the total hydrogen content may be a 
function of the initial orbital period with larger period producing a larger the hydrogen 
content. A more detailed study of this point is left for future work. 

Note that the structure of the donor star at the moment of the explosion is consistent 
with an extended SN lib but with very little H mass (< 0.1 Mq). 
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We analyzed the effect of the secondary star on the observed HST pre-explosion 
photometry. For all the values of at the moment of the explosion of the donor, the 
secondary star is still near the ZAMS. This is a direct consequence of our assumption that 
the object has a mass appreciably lower than that of the donor. The effective temperature of 
the companion is far higher than that of the donor with values within of 22 to 40 thousand 
Kelvin. Thus, the largest contribution to the flux of the system from the secondary is in the 
bluest observed band, F336W, producing a marginal detection of 0.6-2cr level depending of 
the value of ^. Unfortunately, the available HST pre-SN observations are not very suitable 
to constrain the properties of the secondary. 

The ultimate proof of the binary nature of SN 2011dh must come from the possible 
detection of a very hot star once the SN light fades enough. This situation would be similar 
to what occurred with SN 1993J but with different properties of the companion. In any 
case, we should remark that detecting the companion star of SN 201 Idh would provide 
valuable information on the efficiency of the mass transfer process and evolution of massive 
CBSs in general. 
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Fig. 1. — The evolution of the components of a close binary system of solar composition 
stars of 16 Mq and 10 Mq on an orbit with an initial period of 125 days assuming fully 
non - conservative mass transfer (/3 = 0.00). Solid black line represents the evolutionary 
track of the donor star. Dots along it indicate the mass of the star during the RLOFs. The 
corresponding labels are in solar mass units. The star ends its evolution with a mass of 
4.034 Mq with an effective temperature and luminosity compatible with the data observed 
for the object at the position of the supernova SN2011dh. Meanwhile the donor star evolves 
from the ZAMS to pre-SN conditions, the companion star suffers a much slower evolution. 
Its evolutionary track is depicted with a solid blue line and ends with a dot representing 
the conditions attained at the moment of the explosion of the primary. For comparison, the 
ZAMS corresponding to objects of the same composition is shown in dashed line. 
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Logio(pJ [g cm 3] 

Fig. 2. — The evolution of the central temperature as a function of the central density. Solid 
line corresponds to the evolution of the donor star for the case of /3 = 0. Black color indicates 
detached condition, while in red we depict the evolution during RLOFs. It is remarkable that 
during the first RLOF, in which three fourths of its initial mass are transferred, corresponds 
almost to a point in this plane while the other episodes occur on appreciable parts of this 
track (especially the last RLOF). For the sake of comparison, we show in dotted line the 
evolution of an isolated star of the same initial mass {16 Mq) neglecting stellar winds. This 
is very similar to the track corresponding to the donor star. Only minor differences are found 
after carbon core exhaustion. For other values of /3(> 0) the corresponding tracks are almost 
indistinguishable from that corresponding to /3 = and are not included. 
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Fig. 3. — Same as Fig. [T]but for the case of /3 = 0.25. Comparing with the results corre- 
sponding to (3 = 0.00, the star ends its evolution with a mass slightly higher (4.118 Mq), 
while its effective temperature and luminosity are compatible with the data observed for the 
object at the position of the supernova SN2011dh. The evolution of the secondary object is 
represented with a solid blue line. Here, labeled dots along each evolutionary track corre- 
spond to given mass values (in solar units) during the RLOF of the primary. Notice that, 
at the moment of the explosion of the primary, the secondary star is much brighter than in 
the case of (3 = 0.00. 



- 31 - 




Fig. 4. — Same as Fig. [T]but for the case of /3 = 0.50. Comparing with the results corre- 
sponding to P = 0.00, the star ends its evolution with a mass slightly higher (4.118 Mq) 
and its effective temperature and luminosity are compatible with the data observed for the 
object at the position of the supernova SN2011dh. The evolution of the secondary object is 
represented with a solid blue line. Dots along each evolutionary track have the same meaning 
as in the previous figures. Notice that, at the moment of the explosion of the primary, the 
secondary star is much brighter than in the case of /3 = 0.25 and brighter than the donor 
star when it was at the ZAMS. 
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Fig. 5. — Same as Fig. [T]but for the case of /3 = 0.75. Comparing with the results corre- 
sponding to P = 0.00, the star ends its evolution with a mass slightly higher (4.073 Mq) 
and its effective temperature and luminosity are compatible with the data observed for the 
object at the position of the supernova SN2011dh. The evolution of the secondary object is 
represented with a solid blue line. Dots along each evolutionary track have the same meaning 
as in the previous figures. Notice that, at the moment of the explosion of the primary, the 
secondary star is brighter than in the case of /3 = 0.50. 
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Fig. 6. — Same as Fig. [T]but for the case of conservative mass transfer (/3 = 1.00). Comparing 
with the resuhs corresponding to /3 = 0.00, the star ends its evolution with a mass shghtly 
lower (4.014 Mq) and its effective temperature and luminosity are compatible with the data 
observed for the object at the position of the supernova SN2011dh. The evolution of the 
secondary object is represented with a solid blue line. Dots along each evolutionary track 
have the same meaning as in the previous figures. Notice that, at the moment of the explosion 
of the primary, the secondary star is brighter than in the case of /3 = 0.75. 
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Fig. 7. — In order to ease its comparison, we show the evolutionary tracks of the donor stars 
for all the values of (5 considered in this paper. Notice that all tracks are very similar each 
other. 
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Fig. 8. — Same as Fig. [7] but for the case of the secondary star. In sharp contrast with the 
case of the donor star, the evolution of the secondary star is strongly dependent on the value 
of p. We show the ZAMS with dashed line on which we have indicated some values of mass 
(in solar units). For the cases of /3 > the final position in the HRD for the accreting star 
is somewhat hotter and overluminous than objects of the same mass on the ZAMS. 
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Fig. 9. — The evolution of the central temperature as a function of the central density for the 
secondary, accreting star. This is the equivalent to Fig. [2] which corresponds to the central 
evolution of the donor star but, in any case we should remark that the scale is completely 
different because the secondary star is still burning hydrogen during its evolution up to the 
explosion of the donor star. 
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Fig. 10. — The mass transfer rate of mass from the donor star as a function of time. In each 
panel we show the cases for the five values of ^ considered in this work. In all cases there 
occur three RLOFs that have a very similar profile. The only significant change is in the 
time spent by the star from the end of the first RLOF to the onset of the second one. These 
differences are related to changes in the orbital semi-axis. 
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Fig. 11. — The temporal evolution of the orbital period for the five values of /3 considered 
in this work. 
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Fig. 12. — The final outermost hydrogen profile of pre-SN models prior to the explosion for 
all values of /3. Notice that these are very similar each other. 
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Fig. 13. — The spectra of the donor and accreting stars (shown with red an cyan hnes 
respectively). Gray hne represents the addition of both spectra which, in turn, should 
represent the observed one. We show the results corresponding to fully-non conservative 
and conservative cases {(3= and 1 respectively). The mean synthetic (black squares) and 
observed fluxes (red squares) in each bandpass are included in the figure. The secondary 
star is significantly bluer than the primary at the time of explosion. Its contribution to the 
total flux is non-negligible only for the bluest observed band, F336W. See § 13.31 for further 
discussion. 
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Fig. 14. — The red part of the evolutionary tracks for solar composition stars of binary 
systems with masses of 16 Mq + 10 Mq, a fixed value of {3 = 0.50 and different values of 
the initial orbital period. The larger the initial period the lower the effective temperature 
attained by the pre-SN. For Pi = 125 d we find the pre-SN object inside the error box allowed 
by photometry previous to the explosion without need of any fine tuning. 
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Fig. 15. — The evolution of binary systems with a mass of the primary of 16 M©, /3 = 0.50 
and different values for the initial mass of the companion. Labels A and D stand for accretor 
and donor stars respectively. For simplicity, we computed the evolution of the donor with 
the 15 Mq companion and assumed it to be the same for the other companion masses. For 
this donor initial mass value, secondary stars with masses up to ~ 15.25 Mq fall close to 
the ZAMS at the moment of the explosion of the donor star and most of its light will be 
emitted on the blue part of the spectrum. However, if the secondary star has an initial mass 
of 15.50 Mq, at the supernova event it will have at much lower effective temperature. This 
is not compatible with the available observations of the progenitor of SN 201 Idh. 



